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Introduction and Definition 



The improper integral defined by 



/ f(x)dx = lim < / f(x)dx> 

J a b^oo [J a J 



exists when the limit exists. When the limit does not exist, the integral is 
conventionally said to not exist, or to diverge. As an example of this case, 
take f(x) = sin(x). Then 

rb rb 

/ f(x) dx — sin(x) dx = cos(a) — cos(6) 

J a J a 

which oscillates as b — > oo, and does not approach a single value. 

With only the assumption that the integral does exist, however, it's easy 
to show (using sin(x + n) = — sin(x)) that its value would have to be cos(a). 
This suggests that it may be possible to define an integral that applies to 
improper integrals which oscillate as the limit on the upper bound is taken. 

Extending the family of integrable functions to include an additional class 
of functions is on its own sufficient motivation for pursuing this development. 
As a practical motivation, consider a situation where one needs to evaluate 



d_ 

dy 



f(x,y) dx 



and would be able to obtain a greatly simplified integrand if the derivative 
could be brought inside the integral, as 



dx 



POO 


df{x,y) 


la 


dy 



but finds that the resulting integral does not exist. 

The sharp termination of the finite integral at the artificial upper bound 
b in the conventional definition is arbitrary. Further, it is that sharp ter- 
mination that is the source of the oscillations when taking the limit in the 
example above. In order to have a single limiting value, the oscillations that 
occur when taking the limit need to be removed. 
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To this end, we propose the following alternate definition: 

poo [ pb pb+c "1 

f f(x)dx = lim < / f(x)dx + / f(x)z(x — b) dx > (1) 

J a b^oo [J a Jb J 

where f(x) is the function to be integrated, and where z(x) is the termination 
function (defined below). The inclusion of the additional term inside the 
limit allows us to rigorously achieve convergence. To the best of the author's 
knowledge, the definition in (0Q) has never been proposed elsewhere. 

Throughout this paper, for clarity, we include the overstruck Z on inte- 
grals using the proposed definition and a termination function, to distinguish 
them from integrals that exist under conventional definitions. We will assume 
the Riemann definition for conventional integrals when necessary (including 
its extension to improper integrals with an infinite limit), but most of the 
development does not depend on the particular definition. We will point out 
where the Riemann definition is assumed. 

When the integrals and the limit in ([T]) exist, the integral of f(x) is said 
to exist with respect to z(x) under our definition, and to have the value of 
the limit. If, for all termination functions for which the limit in (OQ) exists, 
the value of that limit is the same, then the integral of f(x) is said to exist, 
with that value for the integral. 

We will examine the case where the upper bound is infinite. The corre- 
sponding case where the lower bound is infinite follows easily, and the case 
where both bounds are infinite is found as the sum of two integrals with one 
infinite bound each, in the obvious manner. We will not consider the case of 
improper integrals with finite bounds. 

A termination function is required to satisfy the following restrictions: 

z(x) is finite for all x, (2) 

z(x) = 1 for x < 0, (3) 
z(x) = for x > c, (4) 
and its derivative is required to satisfy 

/oo pc 
z\x) dx = / z'(x) dx = — 1. (5) 
-oo JO 

A termination function is any function satisfying the conditions given in 
(J2]) through A termination function for f(x) is a termination function 
for which the limit in ([T]) exists. 

Implicit in the requirements of ([2]) through (J5j), is that z(x) is not a 
function of b. If z(x) were allowed to vary with b in an arbitrary fashion, 
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for some integrals (where conventional integration fails), one could carefully 
choose z(x, b) to cause the limit in (PQ) to be any desired value. Not all forms 
of variation of termination functions WRT b will necessarily be a problem, 
and some restricted forms of variation WRT b may allow us to expand the 
domain of integrable functions. 

The condition is likewise perhaps not necessary, but makes some of the 
ensuing development simpler. The condition (j3J) also is not strictly necessary. 
A less restrictive condition, z(x) — > as x — ► oo, coupled with the integral 
of z'(x) being absolutely convergent, is sufficient for the developments here. 
For this paper, however, we will proceed using the conditions (J2J) through 
(jSJ), and require termination functions have no variation WRT b. Expanding 
the allowable set of termination functions is a subject for future research. 

The definition in ([T]) can be manipulated into a more useful form. Begin- 
ning with (pQ), 



J. f(x) dx = hm ^ F(x) 



dx + J f(x + b)z(x) dx > (6) 



where 

F(x) = I f(x') dx'. (7) 



Using integration by parts, 

c 

dx — F(x + b)z'(x) dx 



c 

c 











i.°°f(x)dx= lim \F(b)-F(a) +F(x + b)z(x 

Ja b^oo 

Now using conditions (Ej) and (jl]), we obtain 

i°°f(x) dx = -F(a) - lim ( 1° F(x + b)z'{x) dx\ (9) 

J a b^oo [Jo J 

where F(a) is constant WRT b, and is pulled outside the limit. We will make 
extensive use of (9). 

Given z(x), z'(x) is found by differentiation. If we are instead given z'(x), 
z(x) can be found as 

z {x) = [-z'{x)\ ® (1 - s{x)) (10) 
where ® denotes convolution, and where s(x) is the unit step function. Thus, 

/oo /*oo 
z'(x')dx'+ / s(x')z'(x — x') dx' 
-oo J — oo 

poo 

= 1 + / z [x — x') dx' 



o 
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r — oo 

= 1 - / z'(x') dx' 

J x 

z(x) = 1+ r z'^)^' en) 



Taking the derivative of both sides of (jllj) yields the identity z'(x) = z'(x). 

Given two termination functions, a third can be created by convolving 
their derivatives to get the derivative of the third. Given termination func- 
tions zi(x) and z 2 (x), we have 

, ( s / x < , s 

z;-m — ^ n (12) 

,l/ (0 i>Cj v ' 

with 

-^(x) dx = —1. (13) 



Define 

= — z[(x) ® <4(x) = ~ I z i( x ~~ x ') z 2( x ') dx'. (14) 
Using ffT2]) , this becomes 



C2 



z'(x) = - / z'Ax - x')z'Jx') dx'. (15) 
Jo 

The integrand in (I15p is non-zero only when x — x', the argument of z[, is 
between and c\. Applying also that z' 2 is non-zero when < x' < c 2 , it can 
be seen that 

v / x < . s 

z (x) = < „ (16) 
v ' \ x > c v ' 

where c = C\ + c 2 . 

Since ( 1T41) is the convolution integral, using ffT3l for ^! and z 2 immediately 
gives us 

z'(x)dx= / z'{x)dx = — 1 (17) 



which satisfies condition ([5]). Using (1161) in ( flTl) . condition (J3J) is satisfied, 
and then using ffTTj) . we see that condition (T4]) is also satisfied. The condition 
in (T5]) is also easily seen to be satisfied. Thus, z{x) found using (fT4l and ffTTj) 
is a termination function. 

Note that in place of ( fT4l) . the equivalent definition 

z'(x) = - / Cl ^i(x')4(x - x') dx' (18) 
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could have been used. We will call z(x) the combined termination function 
of Zi(x) and z 2 (x), and will call Z\(x) and 22(2;) components of the combined 
termination function z(x). 

Using f fl~4"l) in ( fill) , we can get the combined termination function as 



z(x) 



z[(x' -t)z' 2 (t) dt 



dx'. 



(19) 



For convenience, we will write this relation between the combined termination 
function and its components using the notation 

z{x) = z l {x)Qz 2 {x). (20) 

Some properties of the alternate definition 

We will now show that the integral with respect to a termination func- 
tion Z\(x) gives the same value as integration with respect to a combined 
termination function having Zi(x) as one if its components, with an arbi- 
trary termination function Z2{x) as its other component. Using ffTTl) . the 
limit term from the RHS of ([HI) becomes 



lim 

b— >oo [Jo 



F(x + b)z[(x) dx 



lim 

b— >oo [Jq 



F(x + b)z[(x) dx 



^2 J dec . 



(21) 



The limiting operation is invariant WRT x', and can be pulled inside of the 
integration over x', giving 



z' 2 (x') 



lim 

b— »oo [Jo 



F(x + b)z[(x) dx 



dx' . 



(22) 



The limit as b approaches 00 is the same as the limit as x' + b approaches 
00 for any finite value of x', so (I22p can be written as 



rc2 r ( rc\ 

— \ z 2 (x) lim < / F(x + x' + b)z[(x) dx 

Jo [_b^oo [Jo 



dx'. 



(23) 



The limit WRT b converges uniformly, so the order of the limit WRT b and 
the integration WRT x' can be safely interchanged (See for example |Spi06| 
for this and subsequent interchanges of limiting operations), giving 



lim 

b— »oo [Jo 



z' 2 (x') / F(x + x + fyz'^x) dx 



dx'}. 



(24) 
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Since x' > 0, and since z[(x) = for x < 0, the lower limit of the integration 
over x can be changed to —x'. Since z[(x) = for x > and c — x' >= c\, 
the upper limit of the integration over x can be changed to c — x' . Thus 




F(x + x' + b)z[(x) dx 



dx' 



lim i / z' 2 (x') 

b^oo [Jo 

lim ( / F(x + b) 



F(x + &)^(x — a;') dx 







dx' 
dx 



lim < / F(x + b)z'(x) dx 

b^oo [Jo 



(25) 



We thus obtain, for z±(x) a termination function of f(x) and a component 
of the combined termination function z(x), 

lim ( [ C1 Fix + b)z'Ax) dx\ = lim ( [° F(x + b)z'(x) dx] . (26) 

6-+00 tJo J b^oo [Jo ) 



Equation flzbj) is crucial result. Given a termination function zi(x) for a func- 
tion f(x), and given any other termination function 22 (x), their combined 
termination function z(x) is also a termination function for f(x), and gives 
the same value. That is, for a function /(x) integrable WRT Zi(x), and with 
z(z) determined by 



J a 



dx 



-F(a) 
-F(a) 



lim if f(x + b)z[(x) dx 

b^oo [Jo 

lim < / f (x + b) z' (x) dx 

b^oo [Jo 



(27) 



Given any two different termination functions Z\(x) and z-i (x) for f(x) 
(i.e. for which the integral ([T]) exists) their combined termination function is 
also a termination function for f(x), and the resulting value of the integration 
is equal to that of z\(x). Since the choice of which termination function is 
zi(x) and which is z<i{x) is arbitrary, a function which is integrable WRT some 
termination function z[x) has the same value for all termination functions. 
Thus, if f(x) is integrable under our definition WRT any termination function 
z(x), we can simply say that f(x) is integrable under our definition, without 
further qualification. That is, our definition gives us a unique answer. 

We immediately see that conventional integration with an infinite upper 
bound is a special case of this definition. Simply take as the termination 
function 

z{x) = {l X ~n (28) 
v 1 x > v 1 
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and then 



poo ( pb pb+c 1 

~L f(x) dx = lim < / f(x) dx + f(x)z(x — b) dx> 

J a b^oo [J a Jb J 

= lim < / f(x) dx > 

b— >oo [J a J 

f(x) dx. 



(29) 



Thus our alternate definition of integration gives the same result as con- 
ventional integration, whenever the conventional integral exists. 

Additionally, we quickly obtain that, given any two functions f(x) and 
g(x) integrable under our definition, any linear combination of them is also 
integrable. Assume 



poo 

ff( 

J a 



x) dx = —F(a) + 



and 



poo 
J a 



dx = —G(a) + 

both exist. Then 

poo 

a tf(x) dx + 

J a 



lim 

b— >oo [Jo 



lim 

b— >oo [Jo 



c f 



F(x + b)z'f(x) dx 
G(x + b)z' g (x) dx 



(30) 
(31) 



poo 

(3 Zg(x) dx 

J a 



a \ -F(a) 



+ P 



-G{a) - 



lim { / Fix + b)z'r(x) dx 
)-+oo [Jo 1 

G[x + b)z' g (x) dx 



lim 

b— »oo [Jo 



(32) 



or 



poo roc 

a Zf(x)dx + (3 lg{x)dx = -[aF(a) + (3G{a)} 

J a J a 

- lim fa r Fix + b)z'Ax) dx + (3 P G(x + b)z'(x) dx] . (33) 

b^oo L Jo JO J 



Using (fT^|) . define 



z'(x) 



Z j (y X X ^ Z>g ^ X ') dx' 



(34) 



then using (|2T)]) twice (once with z'^x) = z[(x), and once with z' g (x) = z[(x)), 
fl33l) becomes 



/•OO /'OO 

a tf(x) dx + (3 Zg(x) dx = — [aF(a) + (3G{a)} 

J a J a 



lim 

b—foo [Jo 



c f +c g 



[aF(x + b)+ (3G{x + b)} z'{x) dx \ (35) 



or 

/*oo /*oo /*oo 

a tf(x) dx + f3 tg(x) dx = ~t [af(x) + f3g{x)\ dx. (36) 

J a J a J a 

Thus, our definition satisfies linearity. 
Examples 

The following shows some example evaluations of integrals which do not 
exist under the conventional definition. If conventional integration were at- 
tempted, example 1 would have constant variation as the upper limit of the 
conventional integral approached infinity. Examples 2 and 3 would have 
unbounded variation. 

Example 1: 

f(x) = sin(o;x) 
F{x)= -cos(ax) 



a 

f i/o n 

z(x) = 



1/2 0<x<n/a 
x > ir/a 

z\x) = -1/2 (5(x) + 5(x - n/a)) 



f 

■J a 



oo 

sin(ax) dx 



cos(aa) 



a 



-l/2 1imM ^ ^ (5(x) + S(x - tt/ck)) dx 

cos(aa) i/oi* f cos { a fy + cos(o;6 + 7r) | 



cos(aa) 



a 



Example 2: 



f(x) = x cos(aa;) 
cos(aa;) xsin(ax) 



cr a 

'3/4 < x < n/a 
1/4 ir/a < x < 2-K/a 
x> 2it/a 

z \x) = -1/4 (5(x) + 25(x - n/a) + 5(x - 2n/a)) 



z(x) = < 
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POO 

Z x cos(ax) dx = 

J a 



cos(aa) asm(aa) 



a 2 a 



j. f-1 /cos(a(a; + b)) ^ (x + b) sin(o;(x + b))\ 
b^oo I 4 Jo \ a 2 a ) 

(5(x) + 25(x - ix /a) + 5{x - 2vr/a)) dx 
cos(o;a) asin(cra) 



a 2 a 



+ 



, , ( cos(ab) bsin(ab) 
1/4 hm { + + 

6-+00 y a z a 

^ /cos(a6 + 7r) (7r/a + 6) sin(o;6 + 7r)\ 
\ a 2 a ) 

cos(afe + 27r) (27r/a + 6) sin(a;6 + 27r) | 



a 2 a 



The trigonometric terms within the limit WRT b cancel, leaving 

J°° , . , cos(aa) asm(aa) 

Z xcos{ax)dx = . 

J a a 2 a 

Example 3: 

f(x) = e^ x sin(a::r) 

e f3x 

F(x) = - ^ (flsm(ax) — acos(ax)) 

1+e ln/ a 0<x<n/a 
x > n/a 



z(x = 



z'(x) = --^^ (e^5(x) + 5(x - n/a)) 

roo g/3« 

1 e^ x sin(o;a;) dx = — (/3 sin(o;a) — a cosfaa)) — 

J a a 2 + p 2 

/a e P(x+b) 

(e^ a 5(x) + S(x - n/a)) dx] 



1 + eC 

Working through the algebra, all the terms which vary with b cancel 
completely, leaving 

roo g/^ a 

Z e^ x sin(ax) dx = — ([3 sin(aa) — a cos(aa)) . 

J a a 2 + p 2 
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When /3 — 0, the result is identical to the result from Example 1. Also, 
when (3 < 0, the result matches the result from conventional integration. 



Interchange of order of integration and differentiation 



We will now examine Leibniz's rule for differentiation under the integral 
sign. Under conventional integration, when certain conditions are met, we 
have that 



d \ f 00 



dx 



roo 


9f(x,y) 


Ja 


dy 



dx. 



(37) 



We will show that this relation also holds when one or both integrals require 
our alternate definition. 
Assume 

joo C rc x (y) "1 

2 f(x, y) dx = -F(a, y) - lim \ I F(x + b,y)z x (x,y) dx } (38) 

J a b^oo |yo J 

exists, with f(x,y) finite and continuous WRT y, over a small neighborhood 
of y. We will also assume that we have a termination function that satisfied 



Z x \X, y) continuous WRT y. 



(39) 



Finally, assume 
df(x,y) 



f 

■J a 



dy 



exists, where 



/ F y (x + b,y)z x (x,y) dx } (40) 
Jo J 

9F(x,y) 



dy 



We want to show that 



d 
dy 



roo 

t f(x,y)dx 

J a 



roo 

J a 



df(x,y) 
dy 



dx. 



(41) 



(42) 



We can always choose a constant value c larger than the maximum of the 
values c x (y) needed in (137|) and c x (y) needed in ( 1381) . For convenience, in the 
following we will simply use c for the integration limit, and will only alter the 
upper limit based on how many termination functions are being combined. 
We will make extensive use of ([26]). 



1 While at first this may seem a severe restriction, we can often use (|Ti|) to obtain a 
termination function that satisfies (|39p from one which doesn't. See Example 5. 
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First, looking at the LHS of 



d 
dy 



poo 

1 f(x,y)dx 

J a 



lim 



frf{x,y + h)dx-frf{x,y)dx 
h 



(43) 



lim — 

h^O h 



lim — 

/i— h 



—F(a, y + h) — lim <^ / F(x + b,y + h)z x (x, y + h) dx \ + 
b^oo [Jo J 

F(a, y) + lim ( / F(x + b, y)z x (x, y) dx 
-F(a,y + h) + F(a iy ) - 



(44) 



3c 



Combining the limits on b, 

~F(a,y)-F(a,y + h) 



lim 

h->0 



h 



l r r 3c 

- lim { / (F(x + b,y + h)-F(x + b,y)) 

h b-*oo [Jo 



lim <^ / F(x + b,y + h)z x (x, y + h) z x (x, y) © z x (x, y)dx \ + 
lim ( / F(x + 6, y)^^, y + h) Q z x (x, y) z x (x, y) dx 



(45) 



(46) 



-F y (a,y) - 



lim 



lim 

fe— >oo I Jo 



3c + 6, |/ + h) - F(x + b, y)) 
h 

z x (x, y + h)Q z x (x, y) © z x (x, y) dx 



(47) 



The limit WRT b converges uniformly, and the limit WRT h exists, so we 
can interchange the order of the limits, 



-F y (a,y) 
lim < lim 



3c 



(F(x + b,y + h)-F{x + b,y)) 



o h 

z x (x, y + h)Q z x (x, y) z x (x, y) dx 



(4f 
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-F y {a,y) - 
lim < / lim 



(F(x + b,y + h)- F(x + b, y)) 
h 

z x (x, y + h)Q z x (x, y) © z x (x, y) 



dx 



(49) 



or finally 

d_ 

dy 



POO 

t f(x,y) 

J a 



dx 



3c 



-F y (a, y) - lim \ F y (x + b, y)z x (x, y) z x (x, y) z x (x, y)dx\ . 



(50) 



Now looking at the RHS of fl4z|) . we have 



r 

J a 



or 



POO 

t 



df(x,y) 
dy 



df{x,y) 
dy 



dx = —F y (x, y) — lim ( / F„(x + b, y)z x (x, y) dx\ (51) 



-F y (x, y) - lim ( / F y (x + 6, y)^(x, ?/) z x (x, y) © ^(x, y) dx 



(52) 



Noting that (150]) and (152!) are identical, we obtain P2|) . Since conventional 
integration with an infinite upper bound is a special case of our alternate 
integration, we also immediately obtain 



d_ 

dy 



/ f(x,y)dx\ = Z 

J a J J a 



df(x,y) 
dy 



dx 



(53) 



when the corresponding conventional and alternative integrations exist. This 
allows derivative operators to be "brought inside" an integral operator in 
cases where it is not allowed under conventional integration. This is an 
important result in terms of justifying the use of this definition of integration. 
We also obtain 



d 
dy 



t f(x,y)dx 

J a 



df(x,y) 
dy 



dx. 



(54) 
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It is an open question whether the existence of the LHS of (1421) implies 
the existence of the RHS. An affirmative answer would be an important 
development. 

Similarly, while the existence of the integral on the RHS of ( )42|) does not 
imply the existence of the integral on the LHS (consider f(x,y) equal to a 
constant), it is perhaps possible that there is always some function g(x) such 
that 



f 

J a 



df{x,y) 
dy 



dx 



dx 



d(f(x,y) +g(x)) 
dy 



t (f(x,y)+g(x))dx 

J a 



(55) 



PCX 

f 

d_ 

dy 

Again, this is an open question. 
Examples involving differentiation 

The following examples show passing differentiation inside the integra- 
tion. In the first example, one of the integrals can be performed convention- 
ally, while in the second, both integrations require our alternate definition. 

Example 4: 

cos(xy) 



fv(x,y) 



f(x,y) 

■j 

d ( cos(xy) 
dy \ x 

F y (x,y) 



sm(xy) 



cos(xy) 

y 



z[x) 



zax 



f 1/2 < x < n/y 
I x > ir/y 

4/2 (5(x) + 5(x-7r/y)) 



where |AS72j 



f(x,y) dx 



cos(xy) 



x 



dx = —Ci(ay) 



Ci{z) 



cos(t) 



dt 



The derivative of Ci(x) is given by |Wei08j . Eq. (6) 

dCi{u) 



du 



COS M 



a 
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so 



d f°° cos(xy) ^ dCi(ay) cos(ay) 



dy J a x dy 
From Example 1, 



f 

J a 



d ( cos(xy)\ J°° . . 

' 1 dx = Z — sm{xy) dx 



so 



a dy \ X J Ta 

cos(ay) 

y 



d f°° cos(xy) j°° d f cos{xy)\ 



dy J a x Ta dy \ x 



and we see that (I55|) is satisfied. 
Example 5: 



f(x,y) = sm(xy) 

- cosjxy) 
F{x,y) = 

y 

f y (x,y) =xcos(xy) 

,_, . . cosfxw) xsin(xw) 
F„(x,y) = — ^ + 



y* x 

1/2 0<x<tt/ ? / 
x > n /y 

z x (x,y) = -1/2 (5(x) + 5(x- it /y))E 

3/4 < x < rr/y 
z(x, y) — \ 1/4 -n/y<x<2-n/y 
x > 2n/y 



2 To satisfy (|39p . we could combine z(x,y) with r(x,y) and s(x,y), denned by 
r x (x,y) -- 



s x (x,y) = 



—y/n < x < ir/y 
else 



-(3 * I/max - 2j/)/7T < X < 7r/(3 * 2y) 

else 



where j/ ma x is the maximum y in the domain of interest. The combination with r(x,y) 
gets us a termination function whose derivative is a rectangular pulse of width 2y/n. 
Combining with s(x,y) then gets us a termination function whose derivative WRT x is 
non-zero over a constant domain of 3 * y max /ft, and continuous WRT y, satisfying ()39p . 
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z x (x, y) = -1/4 (5(x) + 25(x - ir/y) + 8(x - 2n/y)) 
From Example 1, 

4°° ■ ( \ j cos(ay) 

Z sm{xy) ax = 

Ja y 



so 



d_ 

dy 



POO 

1 sm(xy) dx 

J a 



From Example 2, 



f 

J a 



'dsm(xy) 
dy 



dx 



d 


cos(ay) 




dy 


y 




cos(ay) 


a sin(ay) 




y 2 


X 



roo 

~t x cos(xy) dx 

J a 

cos(m/) asm(ay) 



comparing, we see that 
d 

dy 

and (H21) is satisfied. 



roo re 

Z sin(xy) dx = ~t 

J a \ J a 



'dsmi^xy) 
dy 



dx, 



Interchange of order of integration 

We will now examine iterated integration of functions of two variables 
with respect to both of the variables, where one of the integrals requires our 
new definition. The development of this section explicitly assumes Riemann 
integration for the other integration, and would need to be modified to remove 
that dependence. 
Assume 

(56) 



roo 

1 f(x,y)dx 

J a 



exists over the domain a < y < (3 for some termination function z(x,y). 
Assume also that w(y) exists over the same domain, and that 



9(x,t) = / f(x,y)w(y) dy 



(57) 



exists over the domain a < t < /3, for a < x < oo. Further, assume that 

poo 

1 g(x,y)dx (5* 

J a 
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exists over the domain a < y < (3 for some termination function z(x,y). It 
may be that (J3B|) and floTl) imply the existence of but that is not proven 
here. We want to show that 



w(y) 
Analogous to ( 



' roo roo rp 

t f(x,y)dx dy= t I w(y)f(x ) y)dy 

J a J J a J a 

0), define 

rx 

F{x,y) = / f(x',y)dx' 



dx. 



and 



G(x,y) = J g{x',y) dx'. 
We will limit ourselves to functions f(x, y) such that 

G(x,t) = f F(x,y)w{y) dy 



(59) 

(60) 
(61) 



(62) 

holds. We will again make extensive use of (1201) and (I26p . and will again 
simply use c for the upper bound of all termination functions. 

Throughout this section we will assume that the integrals WRT y are 
Riemann integrable. This restriction is due to the proof used below, but 
does not appear to be a fundamental limitation of our proposed integral 
definition. A proof without this restriction is desirable, but unknown at this 
time. 

The RHS of (El) becomes 



roo roc roc 

Z [g(x, (3) — g(x, a)] dx = Z g(x,(3)dx—l g(x,a)dx (63) 

J a J a J a 

= G(a,a)+ lim \ / G(x + b, a)z x (x, a) dx > 
b-^oo [Jo J 

G(a, (3) + lim ( [° G(x + b, (3)z x (x, (3) dx\\ (64) 



G(a, a) + lim < / G(x + 6, a)z x (x, a) (g> z x (x, (3) dx > 

b^oo [Jq J 



2c 



G(a } (3) + lim < / G(x + b, (3)z x (x, a) ® z x (x, (3) dx 

b-*oo [Jo 



(65) 



[G(a,a) - G(a,/3)] + 

Mm \ / [G(x + b, a) — G(x + b, j3)\ z x (x, a) ® z x {x, /?) dx \ (66) 
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[G{a,a)-G{a,(3)} + 

lim { / [G(x + b,a) - G(x + 6,, 



M 



:, ( ./•. n J ® z x (x,P) ® n {2^(2;, Tj)} da? f (67) 



where the product symbol n is taken as multiple convolutions, using (1261) M 
times. Equation (1571) holds for any arbitrary M > and set {r^} subject to 
a < Tj < p. Using ( l62l) . ( IBTj) becomes 



[G(a, a) — G(a, (3)\ + 



2c 



- / F(x + b,y)w(y)dy 



lim 

M I 

^(x, a) <g) ^ x (x, /3) ® JJ {^(x, r,)} cfe > . (68) 
j=i J 

We now replace the integral using the standard Riemann integral limit 
definition. That is, we will require that the interval a to j3 be broken into N 
intervals of length Ay iy that each yi fall within the corresponding interval, 
and that the maximum length Ayi approach zero as N approaches infinity, 

= [G(a,a)-G(a,P)]- 



lim 



2c 



lim 



b— >oo Jo N—*oo 



N 



J2F(x + b, y i )w{y i )Ay i 



Li=l 



M 



:,.(•'•■ n)v :,(x,/3) ® ]\ {z x (x, Tj)}dx \ . (69) 
j'=i 



For every iV and set of y,, we can choose M = N and Ti — yi. Thus, for 
every iV and i in (l6"9"j) . we have that 



lim ( / F(x + 6, yi)w(yi) Ay { 

)->oo I Jo 



N 



z x (x, a) <g> z x (x, (3) <g> JJ {^(x, %)} dx }• (70) 



exists. 



We can then interchange the limit over % and the integration in (I69p to 



get 



[G(a, a) — G(a, 
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lim < lim 

b— >oo N—*oo 



N ,2c 



A ' rZC 

/ F(x + b,yi)w{yi) Ay { 



N 

z x (x, a) <g> z x (x, (3) <g> Yl i z x( x > Vj)} dx 

3=1 



(71) 



get 



Using the existence of f!70p again, we can interchange the two limits to 
[G(a,a)-G(a,P)] - 

N ,2c 



lim 



lim \ Y. [ F(x + b, yi)w(yi) Ay { 

>-oo {. =1 Jo 



N 

z x (x,a) <g> z x (x 7 (3) <g> Yl {zx(x,Vj)}dx 



3=1 



(72) 



Using (ED, the LHS of ([59]) becomes 



w{y) 



-F(a, y) - lim ( / F(x + b, y)z x (x, y) dx\ 

b^oo {Jo ) 



w{y)F{a,y) dy 



dy 



lim / w(y)F(x + b,y)z x (x,y)dx 

b^oo Uo 



[G(a, a) — G(a, 



lim 

b—foo (,Jq 



3c 



w(y)F(x + b, y)z x (x, y) ® z x (x, a) (g z x (x. 



13) dx\ 



dy 

dy. 
(73) 



Replacing the integration with the Riemann integral limit definition, 



[G{a, a)-G{a,P)] 



lim 

N^oo 



N 

E 

i=l 



lim 

b— »oo [Jo 



3c 



w(Vi)F(x + b, yi) 
z x {x, yi) ® z x (x, a) <g> ^(x, /3) cfe ^ 



(74) 



Combining the other N — 1 additional termination functions with each i 
using (1261) . we have 



[G(a,a) - G(a,f3)} 



lim 

N— >oo 



' AT 

E 

Li=l 



lim 

b— >oo t_Jo 



3c 



w(yi)F(x + b,yi) 



N 



z x (x, a) ® ^(x, /3) ® H {^(x, cix [> Ay* 
i=i 



(75) 
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or finally 



[G{a,a)-G{a,P)) 



— lim 

TV^oo 



lim ^ V / w(y,i)F(x + b, y { 



Z x yX y OL J 



N 

z x {x,/3) ® Y[ {z x (x,yj)} dx Ayi 



• (76) 



Equations ( 176]) and ( 1721) are the same, so ( 1591) holds. 
Since conventional integration is a special case of this definition, we im- 
mediately get 



and 



rP 

/ w(y) 

J a 



w(y) 



poo 

1 f(x,y)dx 

J a 



f(x,y) dx 



dy 



dy 



roc 
fa 



w{y)f{x,y)dy 



w{y)f{x,y)dy 



dx 



dx 



(77) 



(78 



when the corresponding conventional and alternative integrations exist. 

The existence of the conventional integral on the LHS of ( 1751) may imply 
that using the alternate definition on the RHS is never necessary, but this 
has not been investigated. 

An interesting question is whether the existence of the LHS of ( lo^j) implies 
the existence of the RHS. Intuition suggests that this is the case, but I have 
no proof. It is clearly not true that the existence of the RHS of ( |59|) implies 
the existence of the LHS. Consider f(x,y) constant, where w(y) integrates 
to zero between a and /3. 



Change of variable of integration 



There is a restriction on manipulating integrations under this definition 
that is not present in the conventional definition. While integration by substi- 
tution can be used in some simple cases to change the variable of integration, 
it can not be used in general. To see this, begin with (pQ), rewritten here with 
u as the variable of integration 

/•do f fl3 r/3+7 ~) 

7 f(u)du= lim <^ / f(u)du+ / f (u)C{u - (3) du) . (79) 

J a (l^oo [J a J & J 

Now let 

u = u(x) x = n _1 (n) (80) 
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subject to the restrictions that u(x) be a monotonic function of x, satisfying 



lim u(x) = oo (81) 

x— »oo 

a = u(a) a = n _1 (o;) (82) 
p = u(b) b = u-\f3) (83) 

and let 

du 

g(x) = fWx))-. (84) 
Substituting, the RHS of (1751) becomes 

f /■«" 1 (/3) du ru- l (P+l) d u \ 

^\L ia) mx)) Tx dx+ L w mx))«u( X )-u(b))-dxj. 

(85) 

Using ( 15T1) and (1831) . we can replace the limit (3 —>■ oo with the limit b — ► oo, 
and using ( 1541) . we get 

/ (yf(x) rfx + / g(x)z(x — b, b) dx > (86) 

J a Jb J 

where 

z(z-6,6) = C(w(x)-w(6))- (87) 

While (15fi]) is equal to (155]) . and looks very similar to (pQ), we see that 
z(-) is an explicit function of b. Variation of termination functions WRT b 
is not allowed under our definition, since this can lead to non-unique values 
for the integrals. To maintain equality in going from (1551) to ( 1561) . an explicit 
variation WRT b is required by ( 1571) . The converse of this is that if a change of 
the variable of integration is performed, and a different termination function 
z(x) that is not a function of b is used in (|86p . the value of the limit, and 
hence the integral, may be different. The precise restrictions on variation of 
z(-) WRT b that must be met is a subject for future research, and will have 
a bearing on when change of variable may be performed. 

In the following example, a change of variable using w-substitution is 
performed, affecting the value of the integral. 

Example 6: Consider a square wave function, with value alternating be- 
tween ±1, 

f(x) = 1 -2*mod(La;J ,2) 
where |_^J is the floor function, and choose 

f [x-[x\] mod(LxJ,2) = 
1 ' \ 1 - [x - [x\] mod ([x\ , 2) = 1 • 
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Using @, 

f°° f(x) dx = -F(0) - lim ( 1° F(x + fe)z'(x) dx} 

JO b->oo [Jo J 

We can select c = 1 and use 

z(x) = l/2 0<x<l 

so 

z \ x ) = -1/2 (6(x) + 5(x + 1)) 



and we then get 



roo 



1 f(x) dx =1/2. 
Jo 



Now, attempt a change of variable. Letting 

x = u + a sin(7r«) 
dx = (1 + cot cos(7ru)) (iti 

the integral becomes 

POO 

f f(u + a sin(7ra)) (1 + an cos(iru)) du. 
Jo 

If we restrict a to 



then 

and we get 



Take 
and 



1 1 

< a < - 

IX IX 



f{u + asin(7ru)) = f{u) 



roo 

f f(u + a sin(7ru)) (1 + arx cos(ttu)) du 
Jo 

roo roo 

Z f(u)du + i an f(u) cos(txu) du. 
Jo Jo 



g{u) = a7if(u) cos(7ra) 



G{u) = arx sin {jx (u — [u\ )) . 
We can use c = 1 and z{u) = 1 — u giving 



z'(u) 



-1 < u < 1 
else 
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to find _ 

f g(u) du = 2a. 
Jo 

We then get 

Z f(u)du + ~L anf(u) cos(ttu) du = 1/2 + 2a ^ Z f(x)dx 
Jo Jo Jo 

and using w-substitution to change the variable of integration can be seen 
not be valid. 

In the restricted case where u is linearly related to x, it can be shown 
that the resulting termination function is not a function of b, so a change of 
the variable of integration may be performed. Assume 



u = r + sx x 

with s > 0. Then 



u — r 



a = r + sa a = (a — r)/s (89) 

p = r + s b b = {ft - r)/s (90) 

7 = sc c = j/s. (91) 



Then (185]) specializes to 

f f(P-r)/s 

lim < / /(r + sx)s dx+ 

K03+7-r)/s) ] 

/(r + sx)C(r + sx — (r + s6))s c?x > (92) 



>(P-r)/s 

or, since 6 — > oo as ft — > oo, 



/ /(r + sx) da; + / f(r + 

J a Jb 



sx)z(x — b) dx > (93) 



where 

z(x - 6) = C(ss - sb). (94) 

For this restricted case, the function z(x) is not a function of 6, and 
satisfies the criteria for a termination function. 

Multiple Integration 



22 



A corresponding definition of integration applicable to multiple integrals 
is the next logical development. The goal of this section is to show that a 
multiple integral version of the 1-D definition is feasible. To that end, a 2-D 
definition analogous to the 1-D expression ([9]) will be given here, but a com- 
plete examination of the consequences of the subject will not be attempted. 
The extension of this definition to three or more dimensions is obvious. 

Begin with a two-dimensional analogue to z'(x) satisfying 

z'(r) = for r external to c (95) 
for some simply closed curve c of finite length, and also satisfying 

z\r)d 2 r= ff z\r)d 2 r = -\ (96) 

JjRc 

where R c is region that is the union of the curve c and its interior. These are 
analogous to (j3J) and (j3J). Note that z'(f) is not some kind of derivative of 
some other function z(r). Rather, the notation was chosen simply to match 
that of the 1-D development. 

Then, for any finite length, simply closed curve 7, define T to be 1 in its 
interior, and in its exterior, 

T(fb) = I 1 f inside ° r 011 7( ^' ^ (97) 
I Of outside 7(7% b) 

The curve j(r,b) will be varied so that, in the limit b — > 00, T(f,b) will 
encompass the entire plane. Here, 6 is a parameter, not necessarily a distance. 
A precise definition of what "encompass the entire plane" means may be "the 
radius r(b) of the largest circle centered at the origin and entirely interior to 
7*(r , b) approaches 00 as b — > 00" . We will also require that 7(7% b) fit entirely 
within some circle with a finite radius R(b) for all b. 
Now define, analogously with ( TTUT) 



w(f, b) = -z'{f) <g> T(f, b). (98) 
It can be seen that, for sufficiently large b, 

w(f, b) = 1 r sufficiently far inside of 7(f, b) 
w(f,b) = r sufficiently far outside of j(r, b). 

Finally, the alternate definition proposed for the two-dimensional integral 
over the entire 2-D plane is 

$°°f(x) d2 r = - lim ^JJf{r)w(f, b) rf 2 f| (99) 
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where the same value must be obtained for all possible curve families 7(r, b). 

This last requirement is necessary, since the integral of f{f) is not abso- 
lutely convergent. Particular (perhaps carefully chosen) curve families can 
have finite limits for a given z'(r), while others converge to a different value, 
and still other curve families do not converge at all. Essentially, we want to 
choose z'(r) to ensure convergence, not j(f, b). 

Using this definition, 2-D integrals over unbounded domains of functions 
which are not absolutely convergent are possible, whereas under the conven- 
tional definition they are not |Tre78] . 

While less development has been achieved for the 2-D integral, and much 
work remains, uniqueness and linearity can be shown. Analogous to (Cbl 
convolution of two termination functions yields a third, 



z'(f) 



z\ (f )z'n(f— r 1 ) d 2 f* ' . 



(100) 



The analogue to (T2"oT) can be shown to hold. Assume is satisfied for 
z[(r), that is 



lim 

b— >oo 



f{r)w\{r, b) d 2 r 



lim 

b— >oo 



T(f-f*,b)z[(F) d 2 f 



d 2 r 



(101) 



exists with the same value for all curve families T(f, b). 

We take c 1; c 2 , and c to be circles of the smallest radii r 1; r 2 , and r 12 
respectively, enclosing the non-zero portions of z[(r), z' 2 (r), and z' l2 (r) re- 
spectively. Equation fllOlH becomes 



lim 

b— >oo 



r(f* — r 1 , b)z[(r*) d r* 

r'<ri 



d 2 f 



r" <T2 



z' 2 (r") d 2 r" lim 



JJf{r) Jf^ r{r-7,b)z[{7)d 2 7 c/ 2 r| 



r"<r2 



z' 2 {r") lim 

b^oo 



r'<n 



T(r — r , 6)^ 1 (r ) d f 



d 2 r 



2^" 



d'f 

(102) 



Equation 



exists for all curve families, with the same value. For each 



f", there is a family T(f— r 1 , b) that corresponds to T(r — r 1 + r", b), so we 
can add f" to the argument of T, giving us 



r 1 <r\ 



T(r-f f + f",b)z[(r') d 2 ? 



d 2 r 



d 2 r" . (103) 
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The limit WRT b is uniform, so we can interchange the order of the limit 
and the integration, giving 



— lim 

b— too 



r" <T2 



T(r-r* + r*',b)z , 1 (r') d 2 ? 

r'<ri 



d 2 r 



d 2 r" 
(104) 



Continuing, 



— lim 

6— >oo 



lim 

b— too 



4(f) 

r"<r 2 



r(r — r +r ,b)z 1 (r)dr 

r'<ri 



r" <T2 



r(r — r , b)z x {r - f ) d 2 f 

r'<r±2 



2^/ 



d Z f 



2^/ 



d z r 



d 2 r 



d 2 r 



lim 

b— too 



lim 

b— >oo 

lim 

b— >oo 



r(r — r", 6) 

r'<ri2 



r"<T2 



d 2 ? 



d 2 f 



T(r — r , &)z (r) d r 

r'<Tl2 



lJJr'<ri2 

f(f)w 12 (r,b) d 2 r\ . 



d 2 r 



Finally, then, 



lim 

b— too 



f(r)wi(r, b) d 2 r 



lim 

b— too 



f{f)w\2{r, b) d 2 f 



(105) 



(106) 



Uniqueness and linearity immediately follow from 
same manner as they followed from 0261) in the 1-D case. 

Conclusion 



in essentially the 



The definition of an improper integral with infinite bounds presented 
here is a more general, and more powerful, alternative to the conventional 
definition. The range of functions which are integrable to infinite limits is 
greatly expanded. 

The proposed definition gives the same results as the conventional defi- 
nition when that applies, and preserves uniqueness and linearity. The new 
definition allows interchange of the order of differentiation and integration 
whenever the two integrals exist under the new definition. Also allowed is 
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interchange of the order of integration of iterated integration, again when the 
integrations exist under this definition. The ability to rigorously interchange 
order of integrations, or order of integration and differentiation, in cases 
where integrals under the conventional definition do not converge, provides 
an added tool for manipulation of complicated integrals. 

The definition presented here is distinct from, but compatible with, the 
Hadamard finite part integral [Zwi92j. It can be thought of as providing 
an intermediate level of convergence, less convergent than integrals which 
are absolutely or conditionally convergent, but more so than integrals which 
require finite part integration for removal of an infinite component. Divergent 
integrands which can be handled using the Hadamard finite part definition 
are also divergent under our definition. Conversely, integrands which require 
our definition are not integrable in the conventional Hadamard sense. The 
two definitions can be used together if necessary, with our definition used for 
the evaluation of the finite part, after the infinite term is removed. 

There are several open questions remaining for this definition. It is un- 
known whether differentiation under the integral sign is always possible, or 
whether bringing another integration inside the improper integral is always 
possible. It is also unknown whether the restriction to Riemann integra- 
bility for interchanging the order of integration can be relaxed. Perhaps 
most importantly, the restrictions on termination functions may be stricter 
than necessary, and could be relaxed, further expanding the range of func- 
tions which are integrable under this definition. A restriction on changing 
the variable of integration, not present under conventional integration, was 
demonstrated, although the precise restrictions that must be satisfied are 
unknown. Answering these questions is desirable, and a subject for future 
research. The results presented here, however, are not dependent on these 
answers. 

A further extension of the alternate definition to two-dimensional (and 
higher) multiple integrals was presented, and shown to satisfy uniqueness and 
linearity. While much more development is clearly needed for the multiple 
integral case, enough progress has been made to show that the basic approach 
of this definition is not limited to one-dimensional integrals. 

If one desires to use a conventional definition of the improper integral 
for infinite limits, the development in this paper still allows more freedom 
to interchange order of integration and differentiation. Since the alternate 
definition is defined in terms of standard integrals and limits, and since inte- 
grals using the alternate definition are equal to those using the conventional 
definition when the latter exists, one need only begin with an integral that 
exists under the conventional definition. Changes of order of integration and 
differentiation which lead to integrands which would be nonconvergent may 
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subsequently be performed, using the alternate definition whenever neces- 
sary. All results obtained using the alternate definition remain equal to that 
of the original conventional definition. 



References 

[AS72] Abramowitz and Stegun. Handbook of Mathematical Functions. 
Dover Publications, Inc., New York, 1972. Eq. (5.2.27). 

[Spi06] Jack Spielberg. Interchange of limits and uniform convergence. 

http://math.asu.edu/~jss/courses/fall06/mat472/limit_interchange.pdf, 
2006. 

[Tre78] W. F. Trench. Advanced Calculus. Harper & Row, New York, 1978. 

[Wei08] Eric Weisstein. Cosine integral. Wolffram Mathworld, 2008. 
http: / /mathworld. wolfram.com/CosineIntegral. html. 

[Zwi92] Daniel Zwillinger. Handbook of Integration. Jones & Bartlett, 
Boston, 1992. 



27 



